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Description 
BACKGROUND 

5 [0001] The present invention relates to SAR interferometry, and more particularly, to hyper-precision SAR interfer- 
ometiy using a duai-antenna multi-pass SAR system and processing method. 

[0002] Previous SAR interferometers were either dual antenna systems, such as are described in L. C. Graham. 
"Synthetic Interferometer Radar for Topographic Mapping', Proc. IEEE, Vol. 64. pp. 763-768, June 1974, achieving 
limited relative accuracy, or were derived from post processing of data from pass-to-pass, as is described by R Li and 
10 R. Goldstein, "Studies of Multi-baseline Spaceborne Interferometric Synthetic Aperture Radars', IEEE Trans, Geocsi- 
ence and Remote Sensing, Vol. 28. pp. 88-97, January 1990. 'achieving limited absolute accuracy due to phase am- 
biguities and lack of relative calibration. 

[0003] The synthesis of dual-antenna and multi-pass interferometry is technically difficult because of problems in 
effectively combining the two types of data. A process of simply averaging data from different passes is not effective 

^5 because the phase functions have imprecisely known differences in scale (interferometer baseline), they have unknown 
systematic error terms (at least affine in range-azimuth), they include possibly severe random noise, and, after all of 
these degradations, are only measurable modulo 2n. This also precludes the superficial application of traditional am- 
biguity resolution techniques, for example the independent application of the Chinese Remainder Theorem in each 
range-azimuth cell. Selecting the phase difference basically from one image pair or another, as a function of phase 

20 variance, also fails for the same reasons. 

SUMIMARY OF THE INVENTION 

[0004] The present Invention is an interferometric SAR system and processing method that combines multi-pass 
25 SAR interferometry with dual-antenna SAR Interferometry to obtain elevation maps with accuracy unobtainable by 
either method alone. A single pass of the dual-antenna system provides coarse elevation maps. High accuracy maps 
are obtained through additional passes, with accuracy determined by the number of passes. The processing method 
employed in the interferometric SAR system combines the acquired data to provide a calibrated, high precision, low 
ambiguity elevation map. using approximate least-squares and maximum-likelihood processing methods. 
30 [0005] The traditional dual-antenna SAR interferometer (used as a component of the present invention) collects two 
complex SAR images from slightly different elevation angles on each pass using two antennas on the same platform, 
providing calibrated maps that have coarse precision but are nearly unambiguous because of the small interferometer 
baseline. The multi-pass method collects two or more complex images using multiple passes of a radar platform with 
a single antenna. Alone, the multi-pass method provides much more precise, but ambiguous and uncalibrated, elevation 
3S maps. The present invention combines the dual-antenna and multi-pass techniques to provkJe unambiguous and highly 
precise maps. 

[0006] The present system and method achieves better accuracy and flexibility than conventional approaches. The 
present system and method combines these two techniques to achieve much better accuracy and may provide accu- 
racy and resolution typically associated with low altitude optical stereo mapping systems together with the all-weather, 

40 day-night mapping capability of radar. 

[0007] The dual-antenna interferometer of the present invention provides the basis for a least-squares fit to locally 
affine phase errors resulting from relative-position uncertainties in the multiple flight paths and thus improves relative 
calibration. The dual-antenna interferometer resolves severe phase ambiguities of the inter-pass interferometer inde- 
pendently at each pixel, and thus handles terrain discontinuities, man-made structures. 

45 [0008] Spatial phase unwrapping resolves moderate ambiguities if they are not completely resolvable by the dual- 
antenna interferometer which yields high accuracy in areas of low-to-moderate relief. The elevation separation between 
multiple passes increases geometrically which improves the accuracy for a fixed data collection cost. The elevation 
map is refined iteratively which reduces the impact of calibration errors and outliers. The phase estinf\ation process is 
bootstrapped from small-separation to large separation image-pairs which reduces computatbn. The calibrated com- 

50 piex image-pairs are processed to create 3-D (volume) images, which helps resolve layover and detect concealed 
objects. 

BRIEF DESCRIPTION OF THE DRAWINGS 

55 [0009] The various features and advantages of the present Invention may be more readily understood with reference 
to the following detailed description taken in conjunction with the accompanying drawings, wherein like reference nu- 
merals designate tike structural elements, and in whteh: 
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Fig. 1 shows an operatbnal scenario using a hyper-precision interferometer system and method In accordance 
with the principles of the present invention; 

Fig. 2 illustrates a processing method employed in the system of Fig. 1; 
Figs. 3a and 3b detail the relative calibration processing step of the method of Fig. 2; 
s Fig. 4 details the coarse elevation processing step of the method of Fig. 2; 

Fig. 5 detail the fine elevation processing step of the method of Fig. 2; 

Figs. SaSd show nnaximum-likelihood elevation estimation employed in the processing method of Fig. 2; and 
Fig. 7 details the volume imaging processing step of the method of Fig. 2. 

10 DETAILED DESCRIPTION 

[0010] Referring to the drawing figures, Fig. 1 shows an operational scenario using a hyper-precision interferometer 
system 1 0 and method 20 in accordance with the principles of the present invention. The hyper-precision interferometer 
system 10 comprises a dual-antenna interferometric SAR system 11. The dual-antenna interferometric SAR system 
IS 1 1 is comprised of a SAR radar having two antennas 1 2, 1 3 spaced apart by a predetermined fixed distance.Bo.i • The 
two antennas 12,13 are used to take successive SAR maps of an image scene 1 5 identified as images A-, and A2, A3 
and A4. and A5 and A^, respectively, pass 1: Aq, A^; pass 2: Ag. A3; ... pass k: A2|(.2. Ag^.i. The image scene 15 includes 
a cliff 16 and a building 17, for example. 

[0011] Coarse-precision SAR maps are generated by each pass of the dual-antenna interferometric SAR system 11 
20 by the image scene 15. Multiple passes of the system 10 by the image scene 15 greatly extends the precision of the 
SAR maps using the processing method 20 of the present invention. 

[0012] The dual-antenna system 10 and method 20 provides for calibration and ambiguity resolution not feasible 
with multiple passes of a single-antenna system. Let ^2k^i denote complex SAR image pairs derived from two 

antennas 12, 13 on the kth pass by the image scene 15, where k = 0, . . K -1 . 

2S [0013] A measurement model employed In the method 20 and system 10 of Fig. 1 will now be described. Complex 
images Am(r^) are measured in SAR slant plane coordinates r^ of the SAR radar 11 obtained from different elevation 
angles. These correspond to the different elevation angles of the two antennas 12. 1 3, respectively. An elevation map 
z(x) in ground coordinates x and a volume image A(x,z) are estimated from the complex images. To obtain the elevation 
map, an intermediate variable, the phase difference between the images, Is estimated. 

30 [001 4] The following model is assumed for complex image-pair product, or interf erogram; 

A„ = a(r)expOcPr„(r)] + n„(r) (1) 

^ where n^ is a complex additive noise field. The phase function contains the following components; locally affine 
systematic phase error that may be different for each m, a reference phase term that is random in r but constant in m, 
a phase noise term that is random in both j and m, and phase signal term that Is a locally-linear function of elevation z. 
[0015] The model for the complex image-pair product (interf erogram) is therefore 

AX = "a(r)l^exp(i<Pm^^} + n„„^ (la) 

where n^, ^ additive noise which may depend on a. The model for the phase difference in terms of the elevatk>n map 
is given by 

45 

9^ As) - »^^) (2) 

SO 

where X is the radar wavelength, B^, „ is the elevation separation between the mth and nth Image either in units of 
distance, referred to loosely as the baseline, 6 is the depression angle, ro is the slant range of a central reference point, 
and n is a locally-affine (linear plus constant) phase difference arising from range, angle, doppler, and other offsets 
in imaging conditions. The constant c^ „ is defined to compensate for the difference between two-way propagation of 
^ the single-antenna interferograms and the one-way propagation of the dual-antenna interferogram, so that 
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(3) 



[0016] It is assumed (but is not a fundamental requirement) that the same dual-antenna separation is achieved on 
each pass, so that 



for k = 0, . . K-1 . Other separations B^^ may be arbitrary. 

[001 7] An elevation-dependent correction In coordinates r^ are made to maximize accuracy, making it conform more 
closely to r^. This is an elevatbn map estimation process that is described below. The synthesis of the multi-pass and 
dual-antenna interferometry techniques is technically difficult due to the difficulty in effectively combining the two types 
of data. A process of simply averaging data from different passes is not effective because the phase functions have 
imprecisely known differences in scale (interferometer baseline), they have unknown systematic error terms (at least 
affine in range-azimuth), they include possibly severe random noise and. in addition to all of these degradations, are 
ambiguous in 2n intervals. This also precludes the superficial application of traditbnal ambiguity resolution techniques, 
for example the independent application of the Chinese Remainder Theorem in each range-azimuth cell. Selecting the 
phase difference locally from one image pair or another, as a function of phase variance, also fails for the same reasons. 
The present processing method 20 overcomes these difficulties. 

[0018] Fig. 2 shows the specifics of the processing method 20 employed in the system of Fig. 1. The processing 
method 20 processes the complex image-pairs derived from the dual-antenna interferometer 1 0. The processing meth- 
od 20 comprises four main processing steps including a relative calibration processing step 21 . a coarse elevation 
processing step 22, a fine elevation processing step 23, and a volume imaging step 24. 

[001 9] The relative calibration processing step 21 includes dual-antenna, low frequency phase error processing (step 
25. shown in detail in Fig. 3a) to remove low frequency phase inconsistencies between dual-channel (same-pass) 
image pairs, inter-pass, low frequency processing (step 26, shown In detail in Fig. 3b) to remove low frequency phase 
inconsistencies between inter-pass image pairs, and inter-pass baseline processing (step 27, shown In detail Fig. 3b) 
to estimate the elevation angle separation between each pass. 

[0020] The coarse elevation step 22 includes noise processing (step 72) that adds data from the dual-antenna inter- 
ferometer to reduce the effects of noise, and spatial phase unwrapping (step 73) to resolve minor ambiguities in the 
coarse elevation estimate. The output of the coarse elevation step 22 provides an elevation map. A detailed discussion 
of the coarse elevation step 22 Is provided below with reference to Fig. 4. 

[0021] The fine elevation step 23 includes Inter-pass residual phase processing (steps 81 , 82) to compute the phase 
update attributable to each pass of the interferometer, average scaled residual processing (step 83) to combine these 
multiple phase updates with the proper relative scaling (step 89) into a single accurate elevation update, update 
processing (step 88) to update the elevation map provided by the coarse elevation step 22, and optional iteration 
processing (loop 90) of the updated elevation map to provkie greater accuracy under adverse conditions. The output 
of the line elevation step 23 provides an enhanced elevation map. The details of the line elevation step 23 are provided 
with reference to Fig. 5. 

[0022] The volume imaging step 24 comprises prediction processing (step 92) to predict the phase and magnitude 
of the dominant surface in the image scene 15, subtraction processing (step 92) to subtract the predicted data from 
actual image data, and residual estimation processing (step 95) to estimate the elevation (frequency) of residual data. 
These processing steps are repeated (loop 97) for each pass of the system 10 across the image scene. The output of 
the volume imaging step 24 is a three dimensional SAR Image. The details of the volume imaging step 24 are provided 
with reference to Fig. 7. 

[0023] Details regarding the above four major processing steps 21-24 will be described with reference to Figs. 3-7. 
With reference to Figs. 3a and 3b, they detail the relative calibration processing step 21 of the method 20 of Fig. 2. 
The relative calibration processing step 21 comprises dual-antenna low-frequency processing 25 which comprises 
calculating the interferogram from the complex SAR Image pairs, Ag^, AgK+v as indfcated in step 31. In step 32, the 
affine phase coefficients are calculated from the interferogram of step 31 using phase fitting techniques such as are 
used in conventional SAR autofocus routines. Interferograms calculated In step 31 are added for k = 0 to k = K-1 , as 
indicated in step 33. In step 34, the affine phase coefficients are calculated from the interferogram of step 33 using the 
same procedures employed In step 32. An affine phase correction function is then computed, as indicated in step 35. 
The phase correction function Is then applied to the complex image, as indicated in step 36. Steps 31, 32. 35 and 36 
are repeated for k = 0, ... K-1 . Optionally, if the affine phase coefficients differ significantly for different values of k, the 
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entire set of process steps 31 -36 ar repeated. 

[0024] The inter-pass, low frequency processing (step 27) and inter-pass baseline processing (step 26) are detailed 
in Fig. 3b. The inter-pass, baseline estimation (step 26) comprises calculating the interferogram (step 51). estimation 
of the phase derivatives in range and azimuth for the complex SAR Image pairs, A21C, Aa^+i, as indicated in step 52. 

s and for the sum of dual-antenna interferograms, A, as indicated in step 53. A linear regression procedur is then 
performed on each set of phase derivatives from the estimation steps 52, 53, as indicated in steps 54, 55. The linear 
terms from each respective linear regression step 54, 55 are averaged, as indicated in step 56. This produces an 
estimated baseline ratio. The estimated baseline ratio is multiplied by the dual-antenna baseline (B^^^) as indicated in 
step 57. Finally the multiplied baseline ratio(B2ic,2k+i). the weight (Wphase). baseline estimates from the GPS/INS. and 

10 the weight (Wgps) are combined to produce a weighted average, which results in a final baseline estimate §2k. a ♦i. 
as indicated in step 58. The inter-pass, low frequency processing step 26 is repeated for k = 0, .... K-2. 
[0025] The inter-pass low frequency processing (step 27) processes the bias terms from each of the linear regression 
steps 54, 55 by computing a linear phase reference function in step 61, and then computing phase bias estimates 
(steps 62, 63) for each of the complex SAR image pairs, P^, Agk+a- P^ase bias estimates from steps 62. 63 are 

IS processed to apply affine phase correction thereto, as indicated in steps 64. 65. The phase-corrected complex SAR 
images. Aq ... Agk+i are processed by the remainder of the process steps 22, 23, 24 of the present method 20. 
[0026] Absolute, position is calibrated using the global positioning system (GPS) and/or extemal ground-control 
points. This is common to all mapping systems and is not discussed in detail herein. However, a number of relative 
calibration problems do arise, for which the present invention provides innovative solutions that are described below. 
20 [0027] The dual-antenna low-frequency phase error processing will now be described in more detail. For the dual- 
antenna image-pairs, A^^, A2k+2' ^^e baseline B2i<^k+i known precisely and any locally affine phase difference Pak^k+i 
is precisely attributable to trends in object elevation. It may be possible, however, for small errors to appear in the affine 
term because of motion errors in the two antennas that do not exactly cancel. Estimating and compensation for this 
effect in processing dual antenna image-pairs, Ag^, Ag^+g, is important because subsequent calibration (step 27) and 
2S elevatk>n estimation procedures (steps 22.^23) depend on it. 

[0028] The locally-affine phase function Po.i ^ ^ (steps 32. 33) to the sum 

2JC-1 

30 Z A2k(E)A^^i(l) (5) 

Conventional approaches to this fitting process may be emptoyed such as conventional SAR autofocus and Fourier 
transform-based image registration techniques. 
35 [00i29] The phase of each dual-antenna conjugate product image Ag^AJij^., is normalised (step 35) so that it produces 
that same fitted ph^se component (step 36). This may be accomplished by estlnrjating PgK^+i using the same method 
used in estimating Pq^, and then performing the multiplication 

Relative phase calibration for each conjugate-product image (interferogram) is performed before fusing multi-pass 
interferometric^data.^The process of applying equations (5) and (6) and the inten^ening fit of Pq.i be repeated if 

^ the difference po.i - P2k.2k+i large. 

[0030] Inter-pass baseline processing (step =s 51 -58) will now be discussed in more detail. The baseline separation 
indicated by GPS data, denoted B„^, is not exact and its error may be significant. This type of error causes a scale 
factor uncertainty In z. With multiple passes and dual-antenna data available for each pass, the GPS estinnates of „ 
are refined by performing a weighted-least-squ^res (WLS) adjustment. The dual-antenna baseline Boj is assumed to 

so be equal for each pass. An estimated baseline B2k,2k+2 between each pass is obtained by a least-squares fit between 
the phase of A and A2kAJk+2' phase derivatives of A and A2kA2k'h2 are calculated in range and azimuth. This 

results in two real valued phase derivative images (one for range, the other for azimuth) for each interferogram. Then 
a linear regression is performed between the phase derivative images derived from A and those derived from each 
^ach regression performs a least-squares fit between two phase derivative images (treating each pixel in- 

ss dependently) resulting in a bias term and a slope temi. (The slope term in range and azimuth should be identical except 
for noise so the two are average before further processing takes place.) This slope term is an estimate of the baseline 
ratb and is multiplied by Pqi to get B2kB2k.^2• following WLS estimate Is then obtained 
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A 

s when Wqp3 is Inversely proportbnal to th variance of Bgi^^i^^g ^pj^ is inversely proportional to the variance of 
Reasonable estimates of these variances are typically available. 
[0031] Inter-pass, low-frequency phase error processing (steps 61-65) removes phase inconsistencies and will now 
be discussed in more detail. For the inter-pass image-pairs, the locally-affine phase difference P„ „ has a significant 
random component that Is indistinguishable from object elevation trends. With this type of error, elevation is only de- 

10 termined relative to some unknown reference plane. However, the additbn of dual-antenna data allows the determi- 
nation of this reference plane, or object trend component. It is then easy to pertorm the phase normalization process 
as given in (6), forcing each conjugate-product to have the correct phase trend 



IS 



Q _ ft m,n 



(8) 



consistent with its baseline separation. Conventional approaches can be used to fit „ to the data. For efficiency, the 

phase derivative regression bias term computed above is utilized. It indicates the difference between the linear phase 
20 of inter-pass interferograms and that of dual-channel interferograms. In addition, it may be useful to iterate the esti- 
mation of p„ and described above, because each is affected by errors in the other. 

[0032] Fig. 4 details the coarse elevation processing step 22 of the method 20 of Fig. 2. The coarse elevation process- 
ing step 22 produces elevation maps from the dual-antenna image-pairs. A21C. ^i^+v In the case of multiple passes, 
the coarse elevation processing step 22 averages the dual-antenna results ^^^^^ to enhance the accuracy of the 

2S coarse-precision map, which Is used later for combining and calibrating the inter-pass image-pairs. The interferograms 
are calculated in step 71 after dual channel low-frequency phase corrections have been made to the complex image 
pairs A^if. A2i(^v Normalized interferograms are then added in step 72. and the summation in (5) is repeated after 
phase normalization per equation (6) so that, except for random noise, the product images add in phase and produce 
A. The phase of A is then wnwrapped in step 73 and the elevation is then computed. This produces the coarse elevation 

30 map z^^03„e(r) Only the principal phase is measurable from a single sample. The resulting ambiguity is very minor for 
the dual-antenna products. A two dimensional least-squares phase unwrapping approach is employed to resolve these 
ambiguities. The unwrapped phase, denoted by $9^ is proportional to z. 

[0033] With reference to Fig. 5, the fine elevation processing step 23 that produces fine elevation estimates and will 
now be discussed in more detail. The coarse elevation estimate is used to initialize the fine elevation estimation process- 
es ing step 23. The maximum-likelihood (ML) estimate includes all images jointly, selecting the z that maximizes the 
conditional distribution P(Ao,...A2k.i/z). The fine elevation processing step 23 comprises calculation of the residual 
phase. One way to select the image pairings (values for (m. n)) is indicated in Fig. 5 where the following 2(K-1 ) inter- 
ferograms are calculated. The interferogram A^\^}^^2 '® ^dded to Ag^+iAgi^+a for K = 0, . . ., K-2. Each such interfero- 
gram pair should be identical except for noise, so they are added to reduce noise. Then the next step is to predict the 
40 phase of each interferogram sum Af„A* given the baseline ratio and the phase $oj of the dual antenna pairs, yielding 



0,1 

45 

[0034] More particutariy. interferograms are calculated for complex images pairs A2K. A2k-i>2 and . A2i(^3 in steps 
81, 82. The interferograms are then added in step 83. The phase reference function is renrK3V6d in step 84 using the 
coarse elevation map that is scaled by the baseline ratio in step 89. This produces an interferogram having residual 
phase. Phase unwrapping of the interferogram is then performed in step 85. The phase unwrapped interferogram is 
so then scaled by the baseline ratb in step 86. The scaled interferogram is then averaged over all passes in step 87. The 
elevation is then updated in step 88 to produce a fine elevation map. 

[0035] As the average residual elevatbns are computed (step 87). they are accumulated to form an average residual 
cpo V ^^^^ niay be a weighted-sum based on an estimated phase variance (computed based on interferogram magni- 
tude). 

55 [0036] The initial phase estimate (elevation map) Is updated by 
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^0.1 ^ ^0.1 + ^0.1 



(10) 



This may also be a weighted summation. 

^ [0037] The fine elevation processing step 23 may be repeated (shown by the optional iteration loop between the 
output of the update elevation step 86 and the scale by baseline ratio step 89) to reduce the effects of erors that impact 
the results nonltnearly. For example, ambiguities in „ not resolved correctly on the first pass can sometimes be 
corrected after obtaining the updated phase estimate $oj of equation (1 0). It can also be performed recursively, instead 
of in a batch mode as described above. 

10 [0038] Selection of the baseline separation will now be discussed. The selection of multiple baselines influences 
map elevation accuracy. One strategy is to increase the baseline separation for each pass so that the amount of 
elevation slope that can be sensed unambiguously increases with each pass. This leads to baseline ratio between 
passes of 



where is the design phase noise and Cq is a design constant. With this type of baseline separation, each pass of 
the sensor bootstraps the ambiguity resolution process for the next pass. The separation of the two antennas should 
be large enough to produce a useful map product on a single pass. 

[0039] Referring to Figs. 6a-6d, they respectively show the marginal likelihood function p(A„,A^z) for three inter- 
ferograms and an ML estimate obtained from the product of marginal likelihoods to illustrate the fine elevation process 

23. This maximum-likelihood formulation leads to an equivalent least-squares approach described below. It selects 
the z that minimizes a norm of the residual phase between a model image and the actual images. An approximate 
solution is described below 

[0040] Figs. 6a-6d show marginal likelihood functions and an ML phase estimate in a single pixel, based on the 
product of marginal likelihoods. Note that 



Fig. 6a shows the dual antenna likelihood function poi (z) = p(AoAi/z). Then A^A* is normalized to remove (p„ „, wherein 
^ Po2(2) = P(AoA2/z). shown in Fig. 6b. Then the phase of the normalized version of A^^A* is unwrapped to obtain the 
residual phase „ wherein po4(z) = p(AoA4/z), shown in Fig. 6c. Then „ '® scaled back into the dual-antenna baseline 
so that it can be averaged with the residuals from other inter-pass pairs, shown in Fig. 6d, to produce a joint likelihood 
• P(Ao...A2K.i/z). 

[0041] The volume Imaging step 24 will now be discussed. Previously, each range/azimuth cell was assumed to have 
^0 a contribution from a point on a single-valued surface, whfch could be reconstructed by estimating surface elevation. 
Surface elevation estimation was reduced to a frequency estimation problem which was solved by estimating phase 
differences. Layover effects from severe terrain, vehicles concealed under trees, and manmade structures can produce 
radar returns with contributions from more than one elevatbn in a particular range/azimuth cell. This poses a more 
general problem of estimating the frequency of a superposition of several sinusoids, I.e. spectral estimation. The ele- 
^ vation aperture may not be sufficiently large for conventional imaging. The nonuniform sampling of the elevation data 
(B^ n 's not constant) is an added complication. 

[0042] Fig. 7 details the volume imaging processing step 24 of the method 20 of Fig. 2. For each of the complex 
images Aq. .... ^+2* ^ P*^^' 0 '® selected for volume processing based on magnitude, as indicated in step 91 . Scatter 
signatures are then computed and subtracted from each image in step 92, using stored elevation and complex ampli- 
so tudes 96. Additional scatterers are then detected In step 93. If all scatters are detected in decision block 94, the complex 
amplitude and elevation are estimated for the scatter in step 95 and are stored (step 96). The procedure is repeated 
for each pixel by way of the feedback loop 97 which causes the procedure to detect the scatters for each pixel. 
[0043] More particularly, Fig. 7 shows an iterative procedure for performing the volume imaging step 24 of Fig. 2. 
The volume imaging step 24 is achieved by estimating a single dominant frequency component (step 95), cancelling 
it (steps 92, 96). and then repeating the process. The details of the volume imaging step 24 are as follows. First, pixels 
are selected for volume imaging according to a signal-to-noise ratio (SNR) or magnitude test (step 91). Those with 
sufficient magnitude for volume imaging are processed further by the following iterative algorithm. 
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[0044] Each select d pixel is tested further to detect the presence of additional scatterers (step 93) at an elevation 
different than indicated by the fine elevation processing steps and any previous iteration of volume innaging. This is 
accomplished by retrieving the complex amplitude and elevation for any scatterers previously detected in that pixel 
(step 96). The signatur s for each scatterer is subtracted from each complex image at that pixel (step 92). The resulting 
s residual is tested. If the sum of the magnitude squared residuals exceeds a threshold then an additional scatterer is 
indicated. This is a type of goodness-of-fit test between the model, comprising all detected scatterers. and the obsen/ed 
complex image values. 

[0045] If no addittonal scatterer is detected (step 94) then the process is terminated for that pixel and the next pixel 
is tested (loop 97). If an additional scatterer is indicated (step 94) then its complex amplitude and elevation are estimated 
10 (step 95) using the coarse and fine elevation estinnation processing described above. The estimated amplitude and 
elevation are stored (step 96) and the iterative processing steps are repeated. 

[0046] Thus there has been described new and improved hyper-precision SAP inter-ferometry that uses a dual- 
antenna multi-piass SAP system and processing method. It is to be understood that the above-described embodiment 
is merely illustrative of some of the many specific embodiments which represent applications of the principles of the 
IS present invention. . 



Claims 

20 1. An interferometric synthetic aperture radar system (10) for producing a highly precise three dimensional SAP 
image, said system (10) comprising: 

a synthetic aperture radar system (1 1 ) having two antennas (1 2, 1 3) separated by a predetermined separation 
distance (Bq ^) which generates data comprising complex image pairs (A^, ^2k+^) representative of an image 

25 scene (1 5); and 

processing means for processing (21-24) the data comprising the complex image pairs (Ag^, A2k+i), for cali- 
brating (21 ) the data comprising the complex image pairs (A^, provide calibrated data, for processing 
(22) the calibrated data to generate a low ambiguity SAP elevation map. for processing (23) the low ambiguity 
SAP elevation map to provide an enhanced SAP elevation map, and for processing (24) the enhanced SAP 

30 elevation map to produce a precise three dimensbnal SAP image. 

2. The system of claim 1 , characterized in that said processing means comprises: 

means for processing (25) the data using low frequency phase error processing to remove low frequency 
35 phase inconsistencies between same-pass image pairs (Agif, A^^^)\ 

means for processing (26) the data using inter-pass, low frequency processing to remove tow frequency phase 
inconsistencies between inter-pass image pairs (A^^, ^2^)'* 

means for processing (27) the data using inter-pass baseline processing to provide a coarse estimate of the 
elevation angle separation between each pass. 

40 

3. The system of claim 2, characterized in that said processing means comprises: 

means for processing the calibrated data to add (72) data from the dual^antenna Interferometer (11 ) to reduce 
the effects of noise; and 

45 means for spatially phase unwrapping (73) the calibrated data to resolve minor ambiguities in the coarse 

elevation estimate. 

4. The system of claim 2 or claim 3, characterized in that said processing means comprises: 

so means for providing inter-pass residual phase processing (61) of the data to compute the phase update at- 

tributable to each pass of the interferometer (11). 

means for providing average scaled residual processing (83) of the data to combine these multiple phase 
updates with desired relative scaling into a single update; and 

means for processing the single update to update (68) the elevation map provided by coarse elevation angle 
55 estimation. 

5. The system of claim 4, characterized in that said processing means comprises: 

means for iteratively processing (90) the updated elevation map to provide an enhanced elevation map having 



8 



EP 0 654 681 B1 



greater accuracy. 

6. The system of any of claims 1 - 5, characterized in that said processing means comprises: 

s means for predicting phase and magnitude of a dominant surface in the image scene to provide (92) predicted 

data; 

means for subtracting the predicted data from the image data to provide (92) residual data; and 
means for estimating (95) the elevation of the residual data. 

10 7. A method (20) of processing synthetic array radar data comprising complex image pairs (A2k' ^^a) generated 
by a dual-antenna interferometric synthetic aperture radar system (10) to produce a precise three dimensional 
SAR image, said method (20) comprising the steps of: 

calibrating (21) the data comprising the complex image pairs (A21C. A2k^0 to provide calibrated data; 
IS generating (22) low ambiguity SAR elevation map from the calibrated data; 

generating (23) enhanced SAR elevation map by processing the low ambiguity SAR elevation nnap; and 
generating (24) a highly precise three dimensional SAR image by processing the enhanced SAR elevation 
map. 

20 8. The method of claim 7, characterized in that said step (21 ) of calibrating the data comprises the steps of: 

processing (25) the data using low frequency phase error processing to remove low frequency phase Incon- 
sistencies between same-pass image pairs (Agj^, A2j<+i); 

processing (26) the data using inter-pass, low frequency processing to remove low frequency phase incon- 
2S sistencies between inter-pass image pairs (A2i(> and 

processing (27) the data using Inter-pass baseline processing to provide a coarse estimate of the etevatbn 
angle separation between each pass. 

9. The method of claim 7. characterized in that said step (22) of generating a tow ambiguity SAR elevation map 
30 comprises the steps of: 

processing (72) the calibrated data to add data from the dual-antenna interferometer (11 ) to reduce the effects 
of noise and 

spatially phase unwrapping (73) the calibrated data to resolve minor ambiguities in the coarse elevation angle 
3S estimation step (27). 

10. The method of claim 7 or claim 8, characterized in that said step (23) of generating an enhanced SAR elevation 
map comprises the steps of: 

40 Inter-pass residual phase processing (81 ) to compute the phase update attributable to each pass of the inter- 

ferometer; 

average scaled residual processing (82) to combine these multiple phase updates with desired relative scaling 
Into a single update; and 

processing (88) the single update to update the elevation map provided by the coarse elevation step (22). 

4S 

11. The method of claim 10, characterized in that said step (23) of generating an enhanced SAR elevation map further 
comprises the step (90) of: 

iteratively processing the updated elevation map to provide greater accuracy under adverse conditions and 
to provide an enhanced elevation map. 

so 

1 2. The method of any of claims 7-11. characterized in that said step (24) of generating a highly precise three dimen- 
sional SAR image comprises the steps of: 

predicting (92) the phase and magnitude of a dominant surface in the image scene (15) to provide predicted 
55 data; 

subtracting (92) the predicted data from the image data; 
estimating (95) the elevation of residual data; and 
repeating (97) the above three steps (92. 95) for each pass. 
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Patentanspruche 

1. Interterometrisches Radar-System (10) mit synthetischer Apertur, zur Erzeugung eines hochgenauen dreldlmen- 
sionalen SAR-Bildes. wobei das System (10) aufweist: 

5 

ein Radar-System (11 ) mit synthetischer Apertur. das zwei durch einen vorbestimmten Abstand (Bq i) vonein- 
ander getrennte Antennen (12. 13) aufweist und Daten mit komplexen Blldpaaren (A^.A2k+i) erzeugt, die fur 
eine Bitdszene (15) reprasentativ sind. und 

Verarbeitungsmlttel zum Verarbeiten (21-24) der die komplexen Bildpaare (Ag^.AgK+i) aufweisenden Daten. 
10 zum Kalibrieren (21 ) der Daten mit den komplexen Bildpaaren (Ag^, Ag^+i ), um kalibrierte Daten bereitzustel- 

len, zum Verarbeiten (22) der kalibrierten Daten, urn eine SAR-Hohenkarte mit einer geringen Mehrdeutigkeit 
zu erzeugen, zum Verarbeiten (23) der SAR-Hohenkarte mit gerlnger Mehrdeutigkeit, um eine verbesserte 
SAR-H6henkarte bereitzustellen, und zum Verarbeiten (24) der verbesserten SAR-Hohenkarte, um ein ge- 
naues dreidimensionales SAR-Bild zu erzeugen. 

IS 

2. System nach Anspruch 1. dadurch gekennzetehnet, daB das Verarbeitungsmlttel aufweist: 

ly^ittel zum Verarbeiten (25) der Daten unter Venwendung einer Niedrigfrequenz-Phasenfehler-Verarbeltung, 
um Niedrlgfrequenz-Phaseninkonslstenzen zwischen BikJpaaren (A2ic, A2k+i) aus demselben Durchgang zu 
20 beseitlgen, 

Mittel zum Verarbeiten (26) der Daten. unter Venwendung einer Interdurchgangs-Niedriglrequenz-Verarbei- 
tung, um Niedrigfrequenz-Phaseninkonsistenzen zwischen BikJpaaren (A^i^ A2K+1) at^s verschiedenen Durch- 
gangen zu beseitlgen. und 

Mittel zum Verarbeiten (27) der Daten unter Venwendung einer Interdurchgangs-Basisllnlen-Verarbeitung, um 
25 eine grobe Schatzung der Hohenwinkeldifferenz zwischen einzelnen Durchgangen bereitzustellen. 

3. System nach Anspruch 2, dadurch gekennzeichnet, daB das Verarbeitungsmittel aufweist: 

Mittel zum Verarbeiten der kalibrierten Daten, um Daten aus dem Doppelantennen-lnterferometer (11) hlnzu- 
30 zuf Qgen (72), um die Effekte von Rauschen zu reduzleren, und 

Mittel zum raumlichen Entfalten der Phase (73) der kalibrierten Daten, um gerlngfuglge Mehrdeutigkeiten in 
der Grobhohenschatzung aufzutosen. 

4. System nach Anspruch 2 Oder 3, dadurch gekennzetehnet, daB das Verarbeitungsmittel aufweist: 

35 > 

Mittel zum Bereitstellen einer Interdurchgangs-Residualphasen-Verarbeitung (81) der Daten, um die Phasen- 
aktualisierung zu berechnen. die jedem Durchgang des Interferometers (11 ) zuzuschreiben ist, 
Mittel zum Bereitstellen einer Mittelwertskalierten-Residualverarbeitung (83) der Daten, um diese mehrfachen 
Phasenaktualisierungen mit gewOnschter relativer Skalierung in einer einzigen Skalierung zu kombinieren, 
40 und 

Mittel zum Verarbeiten der einzigen Aktualisierung, um die von der Grobhohenwlnkelschateung gelleferte H6- 
henkarte zu aktualisieren (88). 

5. System nach Anspruch 4, dadurch gekennzeichnet, daB das Verarbeitungsmittel aufweist: 

45 Mittel zum iterativen Verarbeiten (90) der aktualislerten Hdhenkarte, um eine verbesserte Hdhenkarte mit 

einer groBeren Genaulgkeit bereitzustellen. 

6. System nach einem der AnsprQche 1 bis 5. dadurch gekennzeichnet, daB das Verarbeitungsmlttel aufweist: 

50 Mittel zum Vorhersagen von Phase und Betrag einer dominanten Oberflache in der Blldszene, um \A3rhersa- 

gedaten bereitzustellen (92), 

Mittel zum Subtrahieren der Vorhersagedaten von den Bilddaten. um Resldualdaten bereitzustellen (92), und 
Mittel zum Schatzen (95) der Hdhe aus den Residualdaten. 

55 7. Verfahren (20) zum Verarbeiten von Daten eInes Radars mit einem synthetischen Sichtfeld, wobel die Daten kom- 
plexe Bildpaare (A2K. Agk+i) aufweisen, die mit einem interferometrischen Doppelantennen-Radar-System (10) 
mit synthetischer Apertur erzeugt wurden, um ein genaues dreidimensionales SAR-Blld zu erzeugen, wobel das 
Verfahren (20) die Schritte aufw ist: 
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Kallbrieren (21) der die komplexen Bildpaare {P^^ Ai2M) autwelsenden Daten. urn kalibrierte Daten berelt- 

zustellen, 

Erzeugen (22) einer SAR-Hohenkarte mit geringer Mehrdeutigkeit aus den kalibrierten Daten, 
Erzeugen (23) einer verbesserten SAR-Hohenkarte durch Verarbelten der SAR-Hohenkarte mit geringer Mehr- 
deutigkeit und 

Erzeugen (24) eines hochgenauen, dreidimensionalen SAR-Bildes durch Verarbeiten der verbesserten SAR- 
Hohenkarte. 

8. Verfahren nach Anspruch 7, dadurch gekennzeichnet, daB der Schritt (21 ) zum Kalibrieren der Daten die Schritte 
aufweist: 

Verarbeiten (25) der Daten unter Venvendung einer Niedrigfrequenz-Phasenfehler-Verarbeitung, um Niedrig- 
frequenz-Phaseninkonsistenzen zwischen Bildpaaren (Ag^, A^\^^^) desselben Durchgangs zu beseitigen, 
Verarbeiten (26) der Daten unter Venwendung einer Interdurchgangs-Niedrigfrequenz-Verarbeitung, um Nied- 
rigfrequenz-Phaseninkonsistenzen zwischen BiWpaaren {fi^, A2,^^i)ausverschiedenen Durchgangen zu be- 
seitigen, und 

Verarbeiten (27) der Daten unter Venwendung einer Interdurchgangs-Basislinienverarbeitung, um eine Grob- 
schatzung der Hohenwinkeldifferenz zwischen jedem Durchgang bereitzustellen. 

9. Verfahren nach Anspruch 7. dadurch gekennzeichnet. daO der Schritt (22) zum Erzeugen einer SAR-Hohenkarte 
mit geringer Mehrdeutigkeit die Schritte aufweist: 

Verarbeiten (72) der kalibrierten Daten, um Daten von dem Doppelantennen-lnterferometer (11) hinzuzuf Ogen, 
um die Effekte von Rauschen zu reduzieren, und 

eine raumliche Phasenentfallung (73) der kalibrierten Daten, um geringfOge Mehrdeutigkeiten in dem Grob- 
hdhenwinkel-Schatzschritt (27) aufzuldsen. 

10. Verfahren nach Anspruch 7 oder 8, dadurch gekennzefchnet, daB der Schritt (23) zur Erzeugung einer verbesserten 
SAR-Hohenkarte die Schritte aufweist: 

Interdurchgangs-Residualphasen-Verarbeitung (81). um eine Phasenaktualisierung zu berechnen, die jedem 
Durchgang des Interferometers zuzuschreiben ist. 

Mittelwertskaliene-Residualverarbeitung (82), um diese mehrfachen Phasenaktualisierungen mit gewQnsch- 
ter relativer Skaiierung zu einer einzigen Aktualislerung zusammenzufassen, und 

Verarbeitung (88) der einzigen Aktualisierung, um die in dem Grobhohenschritt (22) erzeugte Hohenkarte zu 
aktualisieren. 

11. Verfahren nach Anspruch 10, dadurch gekennzeichnet, daB der Schritt (23) zum Erzeugen einer verbesserten 
SAR-Hohenkarte ferner den Schritt (90) aufweist: 

Iteratives Verarbeiten der aktualisierten Hohenkarte, um eine groBere Genauigkeit unter ungQnstigen Be- 
dingungen bereitzustellen und um eine verbesserte Hohenkarte bereitzustellen. 

12. Verfahren nach einem der AnsprOche 7 bis 11 , dadurch gekennzeichnet, daB der Schritt (24) zum Erzeugen eines 
hochgenauen. dreidimensionalen SAR-Bildes die Schritte aufweist: 

Vorhersagen (92) von Phase und Betrag einer dominanten Oberflache in der Bildszene (15), um Vbrhersage- 
daten zu erzeugen, 

Subtrahieren (92) der Vorhersagedaten von den Bilddaten, 
Schatzen (95) der Hdhe der Residualdaten und 

Wiederholen (97) der oben genannten dret Schritte (92, 95) fur jeden Durchgang. 
Revendications 

1. Systfeme de radar k antenna synth6tlque interf6rom6trique (10) pour produire une image RAAS tridimensionnelle 
de haute precision, ce systdme (10) comprenant : 

un systfeme de radar k antenne synth6tique (11) ayant deux antennes (1 2. 1 3) s^pardes par une distance de 
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s6paratfon pr6d6termin§e (B^ qui g6n6re des donn6es comprenant des paires d'images complexes (Aaip 
Agk+i) repr6sentatives d'une sc6ne d'image (15); at 

un moyen de traitement pour traiter (21-24) les donndes comprenant les paires d'images complexes (/^^, 
Aak+i). pour 6tatonner (21) les donndes comprenant les paires d'images complexes (A^, Agit+i) afin de pro- 
duire des donndes 6talonn§es, pour traiter (22) les donndes dtalonndes de fagon k g6n6rer une representation 
de site RAAS k faible ambiguTt6. pour traiter (23) la representation de site RAAS k faible ambiguTtd pour 
produire une representation de site RAAS am6lioree, et pour traiter (24) la representation de site RAAS ame- 
lioree pour produire une image RAAS tridimensionnelle precise. 

2. Systeme selon la revendication 1 , caracterise en ce que le moyen de traitement comprend : 

un moyen pour traiter (25) les donnees en utilisant un traitement d'erreur de phase de frequence basse pour 
eiiminer des incoherences de phase de frequence basse entre des paires d'images de la meme passe {A^^ 
A2k+i) : 

un moyen pour traiter (26) les donnees en utilisant un traitement de frequence basse Inter-passe, pour eiiminer 
des incoherences de phase de frequence basse entre des paires d'images inter-passes (Aa^, Ajk+i)! 
des moyens pour traiter (27) les donnees en utilisant un traitement de ligne de base inter-passe pour produire 
une estimation grossiere de la separation d'angle de site entre chaque passe. 

3. Systems selon la revendication 2, caracterise en ce que le moyen de traitement comprend : 

un moyen pour traiter les donnees etalonnees de fa^on k additionner (72) des donnees provenant de I'inter- 
ferometre k deux antennas (11) pour reduire les effets de bruit; et 

un moyen pour effectuer un deroulement de phase spatial (73) des donnees etalonnees, pour resoudre des 
amblguTtes mineures dans restimation de site grossiere. 

4. Systems selon la revendication 2 ou la revendication 3, caracterise en ce que le moyen de traitement comprend : 

un moyen pour effectuer un traitement de phase de residu inter-passe (81) des donnees, pour calculer I'ac- 
tualisatton de phase attribuable k chaque passe de I'interferometre (11), 

un moyen pour effectuer un traitement de residu pondere moyen (83) des donnees, pour combiner en une 
actualisation unique ces multiples actualisations de phases avec une ponderation relative desires; et 
un moyen pour traiter I'actualisation unique de fagon k actualiser (88) la representation de site qui est produite 
par {'estimation d'angle de site grossiere. 

5. Systems selon la revendication 4, caracterise en ce que le moyen de traitement comprend : 

un moyen pour traiter (90) de fa9on Iterative la representation de site actualisee, pour produire une repre- 
sentation de site ameiioree ayant une plus grande exactitude. 

6. Systems selon I'une quelconque des revendications 1 - 5, caracterise en ce que le moyen de traitement comprend : 

un nnoyen pour predire la phase et I'amplitude d'une surface dominante dans la scene d'image, pour produire 

(92) des donnees pr6dites; 

un moyen pour soustraire les donn6es predites des donnees d'image, pour produire (92) des donnees de 
residu; et 

un moyen pour estimer (95) le site des donnees de residu. 

7. Precede (20) de traitement de donnees de radar k antenne synthetique, comprenant des paires d'images com- 
plexes (Agk. A2k+i) qui sont produites par un systems de radar k antennes synth6tiques interferometrique k deux 
antennes (10). pour produire une image RAAS tridimensionnelle precise, ce procede (20) comprenant les etapes 
suivantes : 

on etalonne (21 ) les donnees comprenant les paires d'images complexes {A^^, Ag^+i ) pour produire des don- 
nees etalonnees; 

on genere (22) une representation de site RAAS k faible ambiguTte k partir des donn6es etalonnees; 

on gendre (23) une representation de site RAAS ameiioree en traitant la representation de site RAAS k faible 

ambiguTte; et 

on gen6re (24) une image RAAS tridimensionnelle de haute precision en traitant la representation de site 
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RAAS amdlior^e. 

8. Proc6d6 selon la revendication 7, caract6ris6 en ce que I'dtape (21) d'^talonnage des donn^es comprend ies 
dtapes survantes : 

5 

on traite (25) Ies donndes en utilisant un traitement d'erreur de phase de frequence basse pour 6limuner des 
incoherences de phase de frequence basse entre paires d'images de la mdme passe (Agi^, A2ic^t) » 
on trafte (26) Ies donnSes en utilisant un traitement de fr^uence basse inter-passe, pour dliminer des inco? 
h6rences de phase de frequence basse .entre paires d'images inter-passes (A210 ^k^^)'^ 
10 on traite (27) Ies donn^es en utilisant un traitement de ligne de base inter-passe pour produire une estimation 

grossldre de la separation d'angle de site entre chaque passe. 

Proc^de selon la revendication 7, caract6ris6 en ce que I'^tape (22) de generation d'une representation de site 
RAAS k faible ambiguTte comprend Ies etapes suivantes: 

on traite (72) Ies donnees etalonnees pour additionner des donnees provenant de rinterferometre k deux 
antennes (11) afin de reduire Ies effets de bruit; et 

on effectue un deroulement de phase spatial (73) des donnees etatonnees. pour resoudre des amblguTtes 
mineures dans retape d'estimation d'angle de site grossidre (27). 

Procede selon la revendication 7 ou la revendication 8, caracterise en ce que retape (23) de generation d*une 
representation de site RAAS ameiioree comprend les etapes suivantes : 

un traitement de phase de residu inter-passe (81 ) pour calculer Tactualisation de phase attribuable k chaque 
2S passe de Tinterferometre; 

un traitement de residu pondere moyen (82) pour combiner en une actualisation unique ces multiples actua- 
lisations de phases, avec une ponderation relative desiree; et 

un traitement (88) de I'actualisation unique pour actualiser la representation de site qui est produite par i'etape 
de site grossiere (22). 
30 . 

11. Procede selon la revendication 10, caracterise en ce que retape (23) de generation d'une representation de site 
RAAS ameiioree comprend en outre retape (90) suivante : 

on traite de fa^on iterative la representation de site actualisee pour obtenir une plus grande exactitude dans 
des conditions detavorables et pour produire une representation de site ameiioree. 

35 

12. Procede selon {'une quelconque des revendlcations 7-11, caracterise en ce que retape (24) de generation d'une 
image RAAS tridlmenslonnelle de haute precision comprend les etapes suivantes : 

on predit (92) la phase et I'amplltude d'une surface dominante dans la scene d'image (15). pour produire des 
40 donnees predltes; 

on soustrait (92) les donnees predltes des donnees d'image; 

on estime (95) le site de donnees de residu; et 

on repete (97) les trois etapes (92, 95) ci-dessus pour chaque passe. 



so 
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